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We verify the operator product expansion (OPE) of deep inelastic scattering (DIS) on the lattice and present 
first results of a non-perturbative calculation of the Wilson coefficients. 
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1. INTRODUCTION 

The calculation of power corrections in DIS re- 
quires not only the evaluation of the matrix ele- 
ments of higher-twist operators, but also the com- 
putation of the Wilson coefficients beyond pertur- 
bation theory Q. In Q we saw already that the 
scale dependence of the matrix elements of the 
twist-2 operators matches the scale dependence of 
the perturbatively calculated Wilson coefficients 
at most at large scales, so that we may expect 
large non-perturbative effects. A lattice calcula- 
tion of the Wilson coefficients would also save us 
to renormalize the lattice operators. 

We consider DIS off a quark target. This is suf- 
ficient for the calculation of the 'leading' Wilson 
coefficients. The OPE reads 

W^{p,q) = (p\Jf,(q)Jv(-q)\p) + seagull 

= ^a i ^(A- 1 g )0„(A" 1 p), (1) 

n 

where \p) is an off-shell quark state of momentum 
p which is struck by a photon of momentum q, 
and A is a scale parameter. The C n<llu (A~ l q) 
are the Wilson coefficients, and O n (A~ 1 p) are the 
forward matrix elements of the operators 

O lll ... lln =i } Y^D^...D^, (2) 

where Y is a matrix in Dirac space. In the follow- 
ing we will distinguish between bare and renor- 

*Talk presented by D. Petters at Lattice 1998. 



, where 
where 



malized quantities. In the former A = a 
a is the lattice spacing. In the latter A = 
H is the renormalization scale. 

The Wilson coefficients are independent of the 
target. Let us denote the nucleon quantities by 
the superscript N. So if we multiply the nucleon 
matrix elements of the operators (0) by the corre- 
sponding Wilson coefficients, we get the nucleon 
structure functions: 

WtfpM =J2 C n, f >v(A- X q)0W(A- l p). (3) 

n 

The r.h.s. of (j|), being independent of A, can be 
written in different ways: 



C n ^u{aq) On (ap) 

= C n ^ v {vT x q) Z n (jjta) 0i N \ap), 



(4) 



where Z n (p a) is the renormalization constant of 
the lattice operator |^|,[|. As a result, the renor- 
malized and bare Wilson coefficients are related 
by 



C n ,iiv{^ 1 q) = Z n 1 (na)Cn ifl „(aq). 



(5) 



The nucleon structure functions are then most 
conveniently computed from the product of bare 
Wilson coefficients and lattice nucleon matrix el- 
ements. A further advantage of this approach is 
that operator mixing with higher-twist operators 
is automatically included. 
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Figure 1. The eigenvalues w n against n, normal- 
ized so that the largest value is one. 

The plan is now to compute the tensor Wy. v 
and the operator matrix elements O n for a quark 
in the Landau gauge, and to extract the Wilson 
coefficients from this information. 

2. THE METHOD 

We consider a system of quark momenta 
Pmt 1 < m < M with the photon momentum 
being fixed. The Wilson coefficients are indepen- 
dent of p m , so the problem is to solve the N x M 
system of equations: 



Figure 2. The residual error R 2 against n. 
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where Pm are the ensemble averaged (ampu- 
tated) matrix elements of the operators @) be- 
tween quark states of momentum p m , W Pm are 
the corresponding elements of the hadronic ten- 
sor, and the C n are the Wilson coefficients we are 
looking for. The label n runs over the various op- 
erators. Note that Lorentz and Dirac indices have 
been omitted. We shall take all operators with up 
to three covariant derivatives into account. Since 
we do not symmetrize over the indices fi%, /i2j • • •, 
this includes higher-twist operators as well. 



We write 



OC = W. 



as 



(7) 



To compute C, we apply a singular value decom- 
position, which is the standard method of solving 
overdetermined equations W . We write 



O = UwV 1 



(8) 



where U is a column-orthonormal M x N matrix, 
uu is a diagonal N x N matrix, w — diag(w„), with 
positive real eigenvalues w n , f < n < N arranged 
in descending order, and V T is the transpose of a 
column-orthonormal N x TV matrix. The solution 
C is obtained by applying 



l/diag(fM l )C/ 2 



(9) 



to W. To find the vector C of smallest length, we 
replace l/w n by zero if w n is dominated by noise. 



3. RESULTS 

The calculations are done with standard Wil- 
son fermions on 24 3 48 lattices at (3 = 6.2 and 
K — 0T489. So far we have looked at two con- 
figurations. We employed 71 different quark mo- 
menta. For the photon momentum we took 



q 



2 = (7r/2a) 2 = 17.4 GeV 2 . 



(10) 
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Figure 3. The Wilson coefficient c x (aq) of the 
operator X = 'ip'yuDjxip with i ^ j against n. 

In Fig. f we show the eigenvalues w n of w. 
We see a sharp decrease between n f» 40 and 
50 by approximately a factor of 30, and the w n 
for larger n are probably dominated by noise. 
In Fig. 2 we have plotted the residual error 
R 2 = \W - OC\ 2 /\W\ 2 that results from setting 
l/w m = 0, Vm > n, as a function of n. The error 
drops sharply until n reaches w 40 and then stays 
constant, saying that if one makes a fit with too 
few parameters one does not get a good fit, and 
if one fits with too many one 'fits to the noise'. 
So we must truncate the system at n « 40 — 50. 

The commonly used Wilson coefficients (bare 
of all kinematical factors) are obtained by 

c n (aq) = C n ^ u (aq)/C t ™° l ,{aq), (11) 

where C%f? v (aq) is the tree value. In Fig. 3 we 
show the Wilson coefficient c x (a q) [|| of the oper- 
ator X — ip'juDjxijj with i j, again for various 
degrees of truncation. We find a plateau in the 
region where the eigenvalues w n drop to zero, so 
that it does not matter where exactly we truncate 
the system. We have chosen n = 42. Beyond the 
plateau the errors on w n get large. 

As a first step we have verified the OPE. In 
Fig. 4 we show the Wilson coefficients c n (a q) of 
the lower dimensional operators. We compare 
our results with the tree-level predictions. We 
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Figure 4. The Wilson coefficients c n (aq) of the 
first 48 operators. The open circles are our 
Monte-Carlo results. The solid line is the tree- 
level result. The horizontal axis labels the oper- 
ators, including the 16 different T matrices. 

find that the expected structure is very well re- 
produced. In particular, the Wilson coefficients 
which should be zero by symmetry or other argu- 
ments are actually found to be zero. 

As a first application of our method we have 
determined the Wilson coefficient of the opera- 
tor X already shown in Fig. 3. This operator 
gives the first (non-trivial) moment of the unpo- 
larized structure function ||, (x). We find (at 
q 2 = 17.4GeV 2 , cf. ©) c x (aq) w 1.2. If we 
divide this number by the corresponding renor- 
malization constant, Z^ IS (fia), as given in [Q, 
we obtain 

cf^-lg)!^^. (12) 

This is to be compared with the 2-loop result of 
C^in 2 = q 2 ) = 1.01. 
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